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We consider a process of parametric down conversion where the input state is a bosonic 
thermofield vacuum. This state leads to a parametric down conversion, generating an 
output of two excited photons. Following a thermofield dynamics scheme, the input 
state, initially in a bosonic thermofield vacuum, and the output states, initially in vacuum 
states, evolve under a Liouville-von Neumann equation. 
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1. Introduction 

m 

Proposed in 1975, by Takahashi and Umezawa , thermofield dynamics (TFD) is an 
operator-algebraic approach to quantum st at ist ic al mechanic s and a real time for- 
malism to finite temperature quantum field theory^EEEMZI xhe basic ingredients 
of TFD are the doubling of freedom degrees in the Hilbert space H which describes 
the physical system and the building of a finite temperature vacuum, thermofield 
vacuum, by means of a Bogoliubov transformation realized in the zero temperature 
vacuum state defined into the Hilbert space % ® ti, also called Liouville spaced 
This procedure is made in such a way that the expectation value of any operator 
from H in the thermofield vacuum coincides with the statistical mean value. 

TFD has been largely applied in the st udy of finite tem perature systems, hig h en- 
ergy phys i cs l 9 | 10 | ll | 12 | 13 | 14|15 | 16 | 17 | 18 | t tics feo | 21 | 22 | 23 | 24 | 25 | con . 

densed matter physics I " I ' I, among 

Here we consider a process where the input state is a bosonic thermofield vac- 
uum. This state will be lead to a parametric down conversion interaction generating 
photon excitations as output states. In general, this process is described by a way 
where the input state suffering the parametric down conversion is completelly anni- 
hilated generating as output correlated photons with number operator expectation 
values equal to 1, where the output states were initially vacuum state. 
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2 Thiago Prudencio 

Indeed, parametric down conversion is a non-linear process in the light of quan- 
tizing electromagnetic field where a single photon incides on a crystal with second 
order non-linearity generating in the interaction two photons with resulting fre- 
quencies equal to the sum of the input frequency. The output frequencies are called 
generally signal and i dler, the frequency of the input photon is called generally 
pump frequency | 29 | 30 | 

In the process that we will describe, our aim is to obtain as output states one- 
photon number excitations, considering the distribution of the photon number op- 
erator. Following a TFD scheme where the system evolves under a Liouville-von 
Neumann equation, the input state, initially in a bosonic thermofield vacuum, and 
the output states, initially in vacuum states, are changed under the parametric down 
conversion leading the output states one-photon excitations whose caracterizations 
is described by the mean expectation value of their number operators. 

In fact , thermofield dynamics have been used for studying para metr ic amplifi- 
cation ^1 anc j degenerate parametric amplification with dissipation SB. This last 
can treat the case where the bosonic thermofield vacuum evolves under dissipation, 
such that part of the hamiltonian is non-hermitean. 

In our propose, a parametric down conversion of a bosonic thermofield vacuum is 
considered using equilibrium states, such that the Liouville-von Neumann operator 
can be used in an hermitean form. 

2. Thermofield vacuum 

Giving an operator A acting on a Hilbert space H generated by Fock states its 
expectation value in a given ensemble is expressed by 



where p is the density operator in the corresponding ensemble. In thermofield dy- 
namics this expectation value is evaluated by means of the definition of a thermofield 
vacuum |0(/3)), where j3 = 1/T is the inverse of temperature T (ks = h= 1), leading 
to the same result as statistical approach, i.e., 



As a consequence, the thermal vaccuum state |0(/3)) is associated to the density 
operator p. For this reason we need to describe it in a Hilbert space larger than the 
Hilbert space T-L generated by the Fock states |n). Then, the thermofield vacuum 
|0(/3)) is not a vector state in the Hilbert space T-L described by the Fock states |n), 
but a state in another enlarged Hilbert space H®H, where Ti is the Hilbert space 
conjugated to H. In fact, in order to describe |0(/3)) as a vector state, we need to 
double the degrees of freedom of the Hilbert space T-L by a formal procedure named 
tilde conjugation^, creating the space 

Consider the system described by a thermal equilibrium density matrix 



(A) = Tr(pA), 



(1) 



(0(/?)|A|0(/3)) = Tr(pA). 



(2) 



p = e 



(3) 
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Parametric down conversion of a bosonic thermofield vacuum 3 

where Z — Tr(e~@ H ) is the partition function and the energy spectrum E n of the 
hamiltonian H is given by H\n) = E n \n). We construct the thermofield vacuum in 
the space H^H, in terms of the Fock state basis \n) and unknown vectors |c„) € H, 
by means of the following expansion 

|0(/3))=5>«>N- (4) 

n 

In order to find the vectors |c n ), we can use the equation @, where the operator A 
does not act on |c„) vectors, 

m\AW)) =J2(^nK)HA\n). (5) 

m,n 

On the other hand, we can write 

n 

Comparing equations ([6]) and ([5]), we find 

(Cm|Cn) — 2 ^ &mn' (^) 

Thus, by defining the Fock states \n) £ T-L a basis product can be given to % ® "H 
and we can write |c„) vectors as 

|c n ) = -^e-^ B -|n). (8) 



Now, we can write the thermofield vacuum in terms of In, n) £ % ® %, 



w)) = ^Y, e ~ PEn,2 \ n ^)- ( 9 ) 



In fact, the term vacuum is only appropriated because we can define annihilation 
and creation thermofield operators, that commute in the bosonic case and anti- 
commute in fermionic case, 

gj = e- ie tfe ie , (10) 
a = e- l0 ae t0 , (11) 
fit = e ~ tG tfe ia , (12) 

6^ = e- <e 6e i<5 . (13) 
and the thermofield vacuum can be annihilated 

0010(0))= 6,810(0)) =0, (14) 

where G is an unitary operator mixing \n) £ W and \n) g T-L by acting on the Hilbert 
space % as a. two-mode squeezing operator, 

G = i9(f3) (fitgt _ a&) ( (15) 
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where 9 — 9((3) is a parameter related to a thermal distribution, a+ and 6+ are 
creation operators acting on spaces % and %, and a and b are annihilation operators 
acting on spaces H and %, respectively. In the case of a bosonic oscillator system, we 
have 9 = tanh _1 (e _ ' 9LJ ' /2 ), related to a Bose-Einsten statistics, and for a fermionic 
oscillator system we have 9 = tan _1 (e _,3w / 2 ), related to a Fermi-Dirac statistics. 
The term e _zG is a Bogoliubov transformation. Since the Bogoliubov transformation 
is canonical, the corresponding commutations for bosons or anticommutations for 
fermions are preserved. 

The operators and bp , given by equations (flQ|) and (p~2]) , excite the thermofield 
vacuum generating excited thermofield states. 

In the Liouville space H^V. given by TFD, the zero temperature vacuum state is 
given by |0, 0) , associated to a density operator at zero temperature po = |0) (0| € %. 
By applying a Bogoliubov transformation e~ lG on this vacuum state |0,0), the 
thermofield vacuum is generated at a finite temperature T ~ /3 , 

|0(/5)> =e^ G |0,6). (16) 

In TFD, tilde conjugation rules realize a mapping between A operators acting 
on |n) and A acting on \h). These rules are summarized by 

(AB) = AB, (17) 
(zA + wB) = z*A + w*B, (18) 
(it) = (i)t, (19) 

(A) = ±A, (20) 
[A,B} ± =0, (21) 

where the operators A and B act only in the Hilbert space spanned by \n), and A 
and B act only in the Hilbert space generated by |n), where z and w are complex 
numbers, z* and w* are their resp ecti ve complex conjugated. In equation (|20[) . + 
is for bosons and — is for fermions In equation (|21j) + means commutation for 
bosons and — is anticommutation for fermions. 

In the case of a fermionic oscillator, the space T-L (£> T-L is generated from the zero 
temperature vacuum and its excitations, 

|0,0) = 1|0,0), (22) 

|1,0) =at|0,0), (23) 

|0,1) =6t| ,0), (24) 

|1, 1) =3+6+10,0). (25) 

By applying the Bogolioubov transformation e~ lG on the vacuum |0,0) we arrive 
at the following fermionic thermofield vacuum 

|O( ( 0)) = -^(|O,O>+e^|l,l)). (26) 
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From the normalization condition |0(/3)), we derive a partition function Z = 1 + 
e —Pu _ Thj s state can also be written as 

|0(/3)) = cos6»|0,6) + sin6>|l, 1). (27) 

where 

6 = tan-^e-^r). (28) 

The equation (|27[) asserts that the fermionic thermofield vacuum is in the plane 
generated by |0,0) and |1, 1). In fact, the action of the Bogolioubov transformation 
on the vacuum excitations (|22|) . (|23ll . (f24| and (f25|) is given by 

_ ie /|l,l)\ = /cos0-sm0\ f\l,l)' 
\\0,0)J \sin6 cos9 J \ \0,0), 

and It follows that the fermionic thermofield vacuum, eq. (|27[) . is in the plane 
generated by |0, 0) and |1,1) and it corresponds to a rotation of 9, relatively to 
|0, 0). On the other hand, the action of the Bogolioubov transfomation on (f2"5|) and 
([2^|l has no effect, being equivalent to an identity operator in the plane generated 
by |1, 0) and |0, 1) (sec figure[T|). We can use this relation to calculate, for example, 



(29) 



IM) 



0(/3)> 




|0,0) 



Fig. 1. (Color online) Fermionic thermofield vacuum in the plane generated by 1 0, 0) and |1,1). 

the mean value of the number operator 

(oto) - (O(0)|fit6|OG8)) = y^j, (30) 

which is the Fermi-Dirac distribution, where we have agreement with the statistical 
result as given in the equation ([2]) . We can also write ^ 

(at a) = sin 2 9. (31) 

Similar calculations could be done to the case of a bosonic oscillator, case where 
the thermofield vacuum is expressed by 

oo 

10(0)) = Vl - e-P" J2 e-**>, n). (32) 
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6 Thiago Prudencio 

In this case, the state is generated in the subspace generated by the states 
|0, 0), |1, 1), |n, n), ... £ H <Ei H, such that a representation of the bosonic ther- 
mofield vacuum in this subspace is not simple (see figure [2]). In this case, the state 



|3,3> 




W)} 



11,1) 



10,0) 



Fig. 2. (Color online) Bosonic thcrmofield vacuum in the subspace generated by 
|0,0),|1, l),...,|n,n>,... eH®H. 



can be described in terms of hiperbolic functions 

10(0)) = (coslT 1 i?) (tanhtf)™ \n, ft) 



where 



i? = tanh _1 (e-^ /2 ). 



(33) 



(34) 



As said previously, in terms of a Bogolioubov transformation, we can write |0(/?)) = 
e _lG |0,0), where we exchange 9 by fl. In this case, mean value of the number 
operator 

1 



a^a) = (0(fi)\tfa\0(P)) 



ef 1 "-!' 

which is the Bose-Einstein distribution. We can also write 



' a) = sinh 2 



(35) 
(36) 



3. Parametric down conversion 

We con side r a parametric down conversion described by the following hamiltonian 

(n= 1)123 

H = a f a + Pb + c + c + k (aPP + a f bcj . (37) 

In terms of a TFD approach a tilde hamiltonian is also constructed 

H = a* a + Pb + Pc+k (aPP + a f bc] . (38) 
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Parametric down conversion of a bosonic thermofield vacuum 7 

The system evolves according an unitary evolution of the Liouville operator L = 
H — H, given here by 

L = LjQa'a — u>oa/a + ui\Wb — ui\Wb 

+ UJ2C ! C — 0J2&C 

+ k (tftfc - a<b f c + abc f - abc^ (39) 

Consider the initial state 

|^(0)) = |0G8))„|0) 6 |0) C (40) 
where the input state is in the bosonic thermofield vacuum 



|0(/3)} a = \/l - e-P»o J2 e-W»°\n, n) a (41) 

n 

and the output states are both initially into the vacuum |0)f, and |0) c (see figure 
[3]), where 

|O) 6 =|0,0} 6 (42) 
|0) c =|0,6) c . (43) 

By defining 



\0)i 



a 




|0) c 

Fig. 3. (Color online) Scheme of the parametric down-conversion of the bosonic thermofield 
vacuum. 



A„ = e" ' VT - e~^o (44) 

we can write 

\MV) =^A„|n,n) o |0,0> 6 |0 ! 0) c (45) 

n 

The system evolves according to a Liouville- von Neumann equation 

(L-i8t)\il>p(t))=0, (46) 
whose the formal solution is given by 

\Mt))=e-* Lt \M0))- (47) 
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8 Thiago Prudencio 

We can also write the Liouville operator as 

2 

L = y^^u^hj - hi) + k (atfc^ - a6 t c t + h . cj (48) 

i=0 

where the mode frequencies Wo, Wi and uj 2 satisfy 

LU = LUl + LU 2 . (49) 

Under a Liouvillean evolution during a time t and a small coupling term A, the 
initial state 

|V>(0)) = |0(/3))o|0>i|0> 2 
evolves until Xt = 1 leading to the following state 

|#> - , * . 2 (-iJ2 X nV^\n- l,n) o ll,0) 6 |l,0) c 

+ *5^A n v^|n,^l) o |0,l) 6 |0,l) c ). (50) 

n 

Initially the thermofield vacuum has a Bose-Einstein distribution given by the mean 
value of the number state Nq = a^ct, while the expectation value in the other number 
states Ni = b^b and N 2 = &c is zero, 

n 

(Ni) = 0, (52) 
(N 2 ) = 0. (53) 

We turn on the interaction time t until the coupling term A and the time satisfy 
the relation Xt = 1. In this situation, the new expected values of photon numbers 
are given by 

<^°> = ^TV £ ( 2A ^ 2 - nX n) > ( 54 ) 
2(iV )o „ 

(Ni) = 1, (55) 
(N 2 ) = 1. (56) 

Then, before the parametric down conversion of the thermofield bosonic vacuum, 
the mean photon number is increased in each sector measured where there was the 
states |0)& and |0) c , corresponding to one photon in each side. On the other hand, 
the system as a whole is in an entangled state. 

The residual distribution on the side a, where there was initially a bosonic 
thermofield vacuum |0(/3)) can be rewritten in a more simplified form 

(N ) = (N ) -^nm^- 1 -. (57) 
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This shows that the thermofield bosonic vacuum lost in photon number in order to 
create excitations of 1-photon states in the sides b and c. 

As an example, if we consider (Nq)q = 2. This corresponds to T = loq (ln(2)) 1 
or (3 = ln(2)/cjo- It follows that the coefficients p n are given by p n — -±=e~% ln ( 2 ) ; 
and the final distributon in this case will be given by 



E 



2 -nln(2) 



n e 



1 



(58) 



Considering the bosonic thermofield vacuum until 2-photon excitations terms, we 
can write 



|0(/?)) a = V * e ~*? (|0,0) a + e-^°|l,l) a 



Vl - e- 3 ^ 
+ e-^ |2,2) o ) 

and the state is now written as 



(59) 



\Mt)) = 



v/2A? + 4A| 
-iA 2 V^|l,2) o |l,6) 6 |l,0) 



(-iAi|0,l) o |l,6) 6 |l,0) c 



+ i 



Ai|l,0) a + A 2 \/2|2,l) a |0,1) 6 |0,1) C 



(60) 
(61) 



If we project the side a in the state |0, 0) a , we have 

a (o,o\Mt)) =o. 

This implies that this state has no parts in the vacuum and then cannot be pro- 
jected into it. The states have evolved such that the projection into the one-photon 
excitations will lead 

— i\i 



a(0,l\Mt)) 



a(l,0\Mt)) = 



a(2,l\Mt)) 



a(l,2\Mt)) = 



^2X( + AX 2 2 
iXi 

^/2Xf + 4X1 



1,0>6|1,0) £ 



(62) 

(63) 
(64) 

V2AFT4AI 1 - — (65) 

These projections means that all the possible measured values in a will lead to states 
of one-photon excited states in b and c, and these states are completelly separable 



V^Xf + 4A 



0,l) b |0,l) c 



|0,l) b |0,l) c 



|i,6) 6 |i,o) c 
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in the space T-L^T-L. This result is better to explain the previous one In fact, 
the projections 

g <n-l,nhte)= y ~ t J^L |l,0)tll,0)c (66) 



and 



a (n,n- 1|^> = , " V . 2 |0,l) 6 |0,i) c (67) 

Consequently, in the approximation considered the parametric down conversion of 
thermofield vacuum will lead to only excitations of one-photon states. Once specified 
the state of a Fock state T-L®T-L in a, the states in b and c are specified as excitations 
of one photon states in "H ® "H and are not entangled in the space HQ'H. It follows 
that a set of measurements of the residual state in the Fock basis of T-L <S> "H will 
determine clearly the excitations of one-photon distributions in the sides a and b. 



4. Conclusion 

We have considered a parametric down conversion where the input state is given by a 
thermofield bosonic vacuum state on a side a, related to a Bose-Einstein distribution. 
Turning on the interaction until time of interaction and the coupling constant satisfy 
At = 1 and taking into account a small coupling term A, we arrive at the result where 
the output states on sides b and c, initially on vacuum states, are excited in one- 
photon states, leading to number operator expectation values of one photon number 
in each side, b and c. 

On the other hand, the thermofield bosonic vacuum state is not totally anni- 
hilated in the process, resulting in a residual state whose profile of the number 
operator expectation value is different from a Bose-Einstein distribution, due to 
losing of photons in the interaction. 

We shown that a specific measurements on the Fock state of the space H^H 
will lead to a specification of one-photon states distributions in th e sid es b and c. 

Such ideas can also be of relevance on the experimental side 1221, where ther- 
mofield or thermal like states are involved. 
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